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Abstract 

We show that supersymmetric "Dark Force" models with gravity mediation are viable. To 

this end, wc analyse a simple supersymmetric hidden sector model that interacts with the 
visible sector via kinetic mixing of a light Abelian gauge boson with the hypercharge. We 
include all induced interactions with the visible sector such as neutralino mass mixing and the 
Higgs portal term. We perform a detailed parameter space scan comparing the produced dark 
matter relic abundance and direct detection cross-sections to current experiments. 
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1 Introduction 



There has been much interest recently in the possibihty that there exists a hidden sector containing 
a dark matter particle coupled to a hidden U(l) gauge boson (a "Dark Force") having a mass of 
the order of a GeV that kinetically mixes with the photon [l]-[4j. Such a scenario could explain 
many astrophysical puzzles, such as the positron excess observed by PAMELA |5j, ATIC [6|, and 
Fermi [T], or the direct detection and annual modulation signal of DAMA [s], CoGeN T [9 10 
and CRESST [ll] (if one ignores the disput ed [T2||T3) contradiction due to XENON 100 [U] and 
CDMS jlS]). Following from the work of 16 -18 , elegantly simple supersymmetric models in 



the latter category were constructed in [19| and further examined in |20]. However, these works 
emphasized that, in order to obtain such a light hidden sector, supersymmetry breaking effects in 
the visible sector would necessarily be dominated by gauge mediation, in order that the masses 
should be acceptably small. Thus it is natural to ask whether confirmation of these observations 
would be in contradiction with gravity mediation; in other words, whether it is also possible to 
have a gravity-mediated spectrum of particles that can yield similar phenomenology. This is also 
linked to the interesting question as to whether these models can be embedded into string theory: 
such hidden sectors appear very naturally there - see e.g. |2T-29 - but the problem of finding 



see e.g. [21 

gauge mediation dominance over gravity mediation is notoriously difficult to achieve in globally 
consistent models. 

Beyond the dark matter motivation it is also useful to ask what hidden sector models of 
this form coming from string theory are allowed or excluded by current observations. This is 
because, even if the hidden sector does not comprise (all) the dark matter, there are a wealth 
of experiments capable of probing Dark Forces over a very wide range of hidden gauge boson 
mass and kinetic mixing values. Kinetic mixing was considered in the context of the heterotic 



string in 124 30-32 . It has been examined in type II strings in 21-23,33-371; in 23,27 both 



masses and mixings were considered, and it was argued that the Dark Forces scenario could be 
accommodated provided that there is additional sequestering. In this work, we shall consider 
hidden sector models with the particle content and similar couplings to those in [Toj, but argue 
that when we have gravity mediation domination these can still give interesting phenomenological 
predictions under certain mild assumptions, without requiring additional sequestering relative 
to the visible sector. Although we will discuss the possible explanation of the signals found by 
DAMA and CoGeNT, these will therefore not be our primary motivation: rather we wish to 
explore how simple supersymmetric hidden dark sectors with a hidden U(l) can be constrained 
by observations. 

The paper is organised as follows. In section [2] we describe the model of a supersymmetric 
dark sector that we shall be examining. This is followed by a summary of constraints upon hidden 
U(l)s with hidden matter charged under them in[3j There we also include the reach of future fixed 
target experiments, and illustrate these with an investigation of a simple toy model. Section [4] then 
contains the meat of the paper: the results of the parameter search over our supersymmetric dark 
sector model. We include additional technical details in the appendix: the hidden sector RGEs in 
appendix [Aj the spectrum of the model in appendix [B] (including the mass mixing matrix with the 
visible neutralino in B.2); a review of kinetic and mass mixing of a massive hidden gauge boson 
with the hypercharge and Z in appendix [Cj and a description of the Goldstone boson mixing in 
appendix [D| In addition, in appendix [E| we discuss the supersymmetry-induced Higgs portal term 
and the mixing of the hidden and MSSM Higgs fields; we believe that although the existence of the 
term has been known in the literature the mixing terms have not been explicitly given elsewhere. 
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Included is a calculation of the induced coupling of the hidden dark matter Majorana fermion to 
nucleons. 



2 Supersymmetric dark sectors 
2.1 Supersymmetric kinetic mixing 

We shall consider models that interact with the visible sector primarily through kinetic mixing of 
a hidden U(l) gauge field with the hypercharge. Hence we have a holomorphic kinetic mixing Xh 
between hypercharge with coupling gy (and gaugino the bino, b) and hidden gauge superfield 
Xa with coupling g^, (and gaugino written as A) appearing in the Lagrangian density 

The physical kinetic mixing in the canonical basis |23[|36| is then given by 

X = 9Y9hMxh)- (2.2) 

We shall assume no matter charged under both hidden and visible gauge groups so this relationship 
is valid at all energy scales. Since we are considering string-inspired models with a "hidden" U(l), 
that is, without matter charged under both the visible and hidden gauge groups, we shall take 
the value of the holomorphic kinetic mixing parameter to be of the order of a loop factor [23] : 

X'. -I^- (2.3) 

Here k is a number that must in principle be derived from the high-energy string model, but 
depends only logarithmically upon the spectrum thereof. We shall either take it to be equal to 



one, or to vary by at most an order of magnitude from unity 123 , 26 , 27 . We thus have 



^ = 9y9hY^; (2.4) 

the most commonly taken value for x is thus of the order of 10~^, but smaller values correspond 
to decreasing the hidden gauge coupling which may be extremely small in the case of hyperweak 
groups |23j|27j|38j. Henceforth we shall always use the physical mixing x- 

As befits a well-studied subject, there are a variety of notations. In addition to using x^ we 
shall also adopt the notation used in (39) : 

X = - sine = -s, 

cos e = Cf = \J\ — x^^ tan e = t^ = ^ . (2.5) 

Vl - X^ 



Note however that this differs to the expressions in [19], which defines x = — e = —te,s 



X r,- - L 



, , — , although there they write e instead of e (we added the tilde to avoid confusion 

with the above). On the other hand, (40-45 defines 5 = — x- However, a crucial novelty in this 



work is our relation (2.4) which we shall see in section 3.3 will lead to qualitatively different results 
for the cross-sections. 
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2.2 Hidden matter fields 



The model that we shah consider is the simplest possible without adding dimensionful supersym- 
metric quantities. There are three chiral superfields S,H^,H- with and charged under 
the hidden U(l) with charges ±1. These appear in a superpotential with dimensionless coupling 

W D XsSH+H^. (2.6) 

This is inspired by D-brane models where the singlet is essentially the adjoint of the gauge group: 
the superpotential above arises due to the N = 2-like structure, and there is no renormalisable 
singlet potential due to this; alternatively, there may be = 2 supersymmetry of the couplings 
at some scale, although we shall not enforce this. Such hidden sectors from string theory were 



considered in e.g. [25,29 , and the above model was studied with gauge mediation in [19] where 
it was termed a "hidden sector NMSSM," although we have set the cubic singlet term in the 
superpotential to zero. There then exists a global U(l) symmetry under which S and H- are 
charged; string theory will not respect this, and we consider that it shall either be broken at 
higher order in the superpotential or through non-perturbative effects - but we shall assure that 
it will play no role in the following. 

Once we include soft supersymmetry-breaking terms we have the approximate potential for 
the hidden sector 

V =\Xs\WSH+\^ + \SH.\^ + \H+H.\^) 

+ §im\'-\H-\'-o' 

+ {XsAsSH+H_ + ^MaAA + c.c.) (2.7) 



where ^ = — = xidY / gh)gY \ cos2f3. The approximation lies in the D-term potential; the 



«2 

'^<,r — xyyy / yhjyr ^ 
full form is found in appendix [Ej 

A crucial difference for the phenomenology of the model once we consider gravity mediation 
is, however, that the gravitino is not the LSP, and therefore the dark matter can consist of stable 
hidden sector particles. We can thus perform a full analysis of the model, including the visible 



sector and its couplings, using micrOMEGAs [46-50 to determine the relic abundance and direct 
detection cross-sections. 



2.3 Symmetry breaking through running 

Just as in the MSSM the top Yukawa coupling can, through running from the GUT scale, induce 
electroweak symmetry breaking, so in the model we are considering the Yukawa coupling A5 can 
induce breaking of the hidden gauge symmetry. By choosing the soft masses and couplings at the 
MSSM GUT scale we can then find models at the low energy scale with hidden gauge symmetry 
breaking. A priori the independent supersymmetric parameters are x^9h,Xs and the soft masses 
mH^ifns, As and M\ (the hidden gaugino mass) which we can choose at the high energy scale 
and run down. 



Via (2.4) we are asserting a relation between x gh- Thus if we take k = 1 we reduce the 



number of free parameters in the model by one. However, as described above we shall in certain 
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plots allow an order of magnitude variation in k; hence although this does not strictly reduce the 
number of parameters in the model, it does rather constrain them with important consequences. 
Finally we shall make one further assumption about the parameters: we shall take mH+ = ttt-H- 
at the high energy scale. This is motivated by the fields H± being a non-chiral pair (note that 
we are taking no explicit Fayet-Iliopoulos term for the hidden U{1) which would introduce a mass 
splitting). Otherwise we shall scan over the remaining parameters to find interesting models. 

The two-loop RGEs for the model are given in appendix lAl By taking ms > rnH± at the high 
energy scale, the RGEs naturally drive the soft masses for rnr^^ to be negative at low energies, 
triggering hidden symmetry breaking]^ The visible sector coupling via the kinetic mixing then 
determines which field (-ff+ or H^) condenses; without loss of generality we take x to be negative 



and thus i/+ condenses. Defining A = y — m\\^g/g\^ we have the conditions for a stable 
minimum with = A/A5, all other expectation values zero: 



<m?L + + m| + 2l\^ 

<{m^_ +m\ + A2)(m| + A^) - \As\^^^. (2.8) 
This is reviewed in appendix [Bj The hidden gauge boson mass is then given by 

my = {^J2gh/\s)^. (2.9) 
We give two examples of the values obtained scanning over ms and 05 = ^ in figure |l 

2.4 Symmetry breaking induced by the visible sector 



The mechanism for hidden gauge symmetry breaking promoted in work such as 19 is via the 
effective Fayet-Iliopoulos term induced in the hidden sector by the kinetic mixing with the visible 
Higgs D-term. In such a case, the mass-squareds m^, m?_ may be positive provided they are small 
enough that A^ > 0. One motivation for this work is that such a case is difficult to justify in the 
case of gravity mediation; it is however not impossible, for example through fine-tuning of the 
parameters or through sequestering. We shall in the section |4] examine this case, still assuming 
that the gravitino is not the LSP. 



2.5 Dark matter candidates 

The model above contains essentially two different dark matter candidates: a Majorana fermion 
and a Dirac onej^ Neglecting the effect of kinetic mixing with the visible neutralino, the fermion 
mass matrix in the basis {X,h-^,h^,s) corresponding to hidden gaugino, hidden higgsinos and 
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^ We ignore the effect on the running of the kinetic mixing, since such terms always enter suppressed by 0{x^ 
Of course, it would be interesting to include all of these effects, but the hidden sector running would then be (of 
course very weakly) dependent upon the visible sector parameters and we leave this to future work. 

^We are ignoring the possibility of scalar dark matter since, although the model as we have written it contains 
stable scalars, we expect the symmetries protecting them to be broken at some higher order in the potential allowing 
them to ultimately decay. 
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Figure 1: Hidden photon mass induced by radiative hidden gauge symmetry breaking, scanned 

over ms and as = In both, mu = As = 100 GeV, ah = 0.0417. Left: Ma = 71 GeV, 
right: M\ = 50 GeV. Ah values given at 10^^ GeV. The black region shows no stable symmetry 
breaking. 



hidden singlino is given by 



f Mx rriY \ 



Mi 





V 




A 
AO/ 



(2.10) 



The Majorana particle is formed from diagonalising the A, states; in the case of a large Mx 
this leads to a see-saw effect and a low mass. We shall refer to this state as "oi", micrOMEGAs 
notation for the lightest odd particle. Clearly there will therefore always be a fermion lighter 
than the hidden gauge boson (to avoid this fate we would need to add a mass for the hidden 
singlino). In order for the Dirac fermion formed from s to be the lightest state we would need 
Xs < V^Oh and for the Majorana mass Mx to be rather small at the high energy scale (this could 
happen for example in a string model where the modulus corresponding to the gauge coupling 
does not obtain an F-term) although it is somewhat suppressed in running down to the low scale. 
Hence the Dirac fermion scenario is not compatible with radiative-breaking models, but presents 
an attractive candidate for the visible sector induced breaking. We shall refer to this state as 



11 ~ 5) 
07 . 



In a complete analysis including the couplings and annihilation cross-sections, it is necessary 



to take the mixing with the visible neutralino into account; this we do in appendix section B.2 
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3 Constraints and discovery potential 



There are already a wealth of constraints on the parameter space of models with dark forces and 
hidden matter that we must apply in our search over models. However, there are also future 
experiments which will have the potential to rule further regions out - or make a discovery. In this 
section we summarise these current and future constraints, and illustrate them by application to 
a toy model. 



3.1 Limits on the hidden photon 



Constraints from electroweak precision tests (EWPT) are used as have been presented in 39 
where the strongest constraint is provided by the mass of the Z for most of the parameter space. 
In the following plots of x vs ra^i this is shown as a long-dashed approximately horizontal blue line 
excluding roughly X ^ 3 x 10^^. Another constraint comes from the muon anomalous magnetic 
moment |51j and is dominant for m^i < 1 GeV: in most of the following plots of x vs m^i this is a 
long-short-dashed brown line at low masses and x > 10~^. There are model dependent constraints 
from BaBar searches |39j that might be the most constraining in the region 0.2 GeV < m^i < 
10 GeV but only apply if the 7' can not decay into hidden sector particles; in most following plots 
of X vs this is a dashed dark purple line at low masses below 10 GeV and x ~ 2 x 10~^. These 
do apply for most of the supersymmetric models we are considering, where the mass of the 7' 
and hidden matter are similar - preventing a decay to the hidden sector. However, if the hidden 
photon can decay to hidden matter, then there is instead a much weaker constraint from the Z 
width; we require 

r(Z chidden) ^ , , 

which for a single hidden Dirac fermion of mass Mx < Mz and unit charge under the hidden 
U(l) corresponds to (see also [52j) 



9l\ , M^x 



4^ < 0.008 



(3.2) 



where cw, sw are the usual cosine and sine of the weak mixing angle respectively, is defined in 
equation ( |C.7 ); for Mx ^ Mz this simplifies to X9h ^ 0.04. Clearly for a small number of hidden 
particles (and gh < 1) this is a weaker constraint than the measurement of the Z mass. 
For m. 



53 



54 



below 1 GeV there are additional constraints from fixed target experiments 
Not only are there limits on the kinetic mixing for very light hidden photons, but excitingly there 
are also dedicated experiments planned (and partly already running) that can further probe this 
parameter space with real discovery potential. The estimated sensitivities of those experiments 
are shown in the plots for the toy model in section 3.3 There are two fixed target experiments 
(MAMI (55] and MESA) in Mainz, three (APEX (56|[57], DarkLight [5| and HPS) at JLab and 
one (HIPS) p] at DESY. 



3.2 Limits from dark matter 

There are further constraints arising on the dark matter particle, its mass and its interactions. 
First of all, the dark matter particle should not have a relic abundance in excess of the one 
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measured by WMAP 60 . This is a very strict limit and translates to a limit on the dark 
matter (DM) annihilation cross section. We compute the dark matter relic abundance using 
micrOMEGAs where we have implemented our model. However, while there is an upper limit on 
the relic abundance, there is no objection to having a dark matter candidate whose abundance is 
lower than the one measured. In this case, it would than only be a part of the total dark matter 
(we shall refer to this as subdominant DM) and the remaining dark matter density would consist 
of other particle(s) such as an axion or axion-like particle whose phenomenology is not the subject 
of this article - we shall simply assume in such cases that the direct detection cross-sections and 
interactions with the hidden sector of the additional dark matter are both negligible. In our plots 
we show parameter points that give an abundance in agreement with the WMAP value in dark 
green and ones where the DM is subdominant in light green. 

Additional constraints apply to the dark matter particle and its scattering cross section on 
nuclei. It is necessary to distinguish spin-dependent (SD) and spin-independent (SI) scattering. 
Dependant on whether the dark matter particle is a Majorana or Dirac fermion it has either 
dominantly SD or SI interactions respectively. Current limits from direct dark matter detection 
experiments are strongest for SI scattering cross sections (~ 10~^^ cm^) while SD cross sections 
only start to be excluded at the 10~^^ cm^ level. 

For the low dark matter masses (~ 10 GeV) we are interested in, the most relevant constraints 
come on the SI side from XENON and CDMS. However due to the signal claims from DAM A 
and CoGeNTj^ there has been a large debate on the reliability of those constraints, especially at 
low dark matter masses close to the energy threshold of the experiments. There are also large 
astrophysical (halo model, dark matter velocity, local dark matter density) and nuclear physics 
uncertainties that should be taken into account. Even though XENON and CDMS claim to rule 
out most of the DAMA and CoGeNT preferred regions, the positive signals remain and there have 
been various studies of how to reconcile those different results. |^ We adapt the analysis of 66 



which made a systematic scan taking into account the various uncertainties. There it is found 
that depending on the halo model some of the CoGeNT and sometimes even DAMA preferred 
region is consistent with the exclusions from XENON and CDMS. For the details of the different 
halo models see [66]; we will mostly use their so-called SMH (Standard Model Halo) and in a few 
cases show the differences that arise when changing for example to a NEW or an Einasto profile. 

We strictly apply the XENONIOO and CDMSSi constraints derived in [66] to the SI scattering 
cross sections and only show points that are not excluded by any of the two experiments. In the 
plots of (Tsi vs rriDM in section [4] (see figures [7, 11 and 12) the CDMS limit is shown as a dashed 



turquoise line while XENON is a blue dash-dotted line (in the plots of section 3.3 the additional 
constraints due to accelerator experiments are shown as grey areas) . For most halo models CDMS 
is more constraining at lower masses than XENON. 

Regarding SD scattering there are several direct detection experiments sensitive to the low 
dark matter masses we are interested in. Until June, Picasso for the lightest and COUPP for 
the slightly larger masses were the most constraining experiments |67| . There is also a quite 
strong limit from Super-K using neutrino fluxes produced by dark matter annihilation in the Sun 



^We have not explicitly included the CRESST signal in our search. One of their two signal regions is roughly 
compatible with both DAMA and CoGeNT signals, although this is still subject to astrophysical uncertainties. 
■^One interesting possibility is to allow isospin-dependent interactions with just the right behaviour to suppress 



the interaction cross-section with xenon nuclei 61-65]. We simply note that, although in the case of hidden Dirac 



fermions the interaction is almost entirely with protons rather than neutrons, in our models this tuning is not 
possible. 
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which however only apphes to dark matter masses above 20 GeV (only neutrino-induced upward 
through-going muons have been used in this analysis which leads to a quite high energy threshold 
and therefore a sensitivity only to larger DM masses) |68| . Very recently a new direct detection 
experiment SIMPLE |;69i| has published a limit on the SD scattering cross section which in the 
low mass range is one order of magnitude stronger than previous experiments (for a critique of 



their limit see 70 and the collaboration's response [71]). 

In the following analysis we use the pre-SIMPLE constraints as strict exclusions and show only 
points consistent with those limits. As there has been criticism of SIMPLE's limit we will not 
apply this universally but rather show how our results change when taking it into account. In the 
SD plots the exclusion lines for the different experiments are as follows: SIMPLE brown short- 
dashed line, Super-K black long-dashed line, PICASSO orange almost solid line, COUPP2011 
turquoise dash-dotted line, COUPP2007 blue solid line, KIMS green long dashed line. 

In the plots the scattering cross section on protons both in the SD and the SI is given. For 
SI scattering the reason for this is that the dark matter couples via the kinetic mixing and thus 
to the charge of the nucleon. This is not the case for SD scattering and both cross sections on 
neutrons and protons are roughly the same, however the one on protons is shown. 

Those constraints on the scattering cross section can strictly only be applied to particles that 
actually constitute the entire dark matter density. If the dark matter is subdominant however, 
the limits on its scattering cross section have to be rescaled accordingly: the local density of a 
dark matter candidate ■0 relates to the local total DM density pdm as their abundances 



and so do the limits that are set by direct detection (DD) experiments. Thus an experimental 
bound on ctdd translates into an actual bound on the scattering cross section of if) as 

<TV. = 0-DD7S • (3.4) 

"DM 

This means that direct detection constraints on the scattering cross section become less potent 
for subdominant DM particlesj^ 



3.3 Application to toy model 

To illustrate the above constraints/future experimental reach, and more importantly provide a 
comparison to the more complete model of section [2] that we shall investigate in section [4j here 
we shall consider a toy model. This is the simplest possible dark sector: a Dirac fermion ^ with 
unit charge only under the (massive) hidden U(l). We shall not include any Higgs sector - the 
U(l) could, after all, naturally have a GeV scale mass via the Stiickelberg mechanism p3l i [27| - so 
we will not consider how the dark matter particle becomes massive. This is essentially the model 



considered in [41-45 except that we shall insist on the relation (2.4); the parameters are the dark 
matter mass mDMj hidden photon mass m^/, kinetic mixing x and the tuning parameter k. 



9 




Figure 2: Annihilation diagrams: s-channel annihilation on the left, resonant at my 
t-channel on the right, accessible and dominant when > Tny. 



2mDM; 



3.3.1 Constraints and future searches 

The DM can annihilate through and/or into hidden photons according to the diagrams shown in 
figure [2} Whereas the left diagram is possible for all DM masses, the right one is kinematically 
only accessible when ttidm ^ nT'-y' ■ The left diagram also leads to a resonant enhancement of the 
annihilation cross section and accordingly to a dip in the relic abundance for my = 2mDM- This 
can been seen in figure [3] where we show the relic abundance for a dark matter mass of 6 GeV 
(two top plots) and 7 GeV (bottom) in function of the kinetic mixing x and the hidden photon 
mass m^/. The grey areas are those excluded by beam dump experiments (the curves on the left 
side of the plot) and muon and electron anomalous magnetic moment (top left corner of the plot) 
and the BaBar search as well as EWPT for large kinetic mixing and hidden photon mass. The 
thin dark green band is the region which gives the correct WMAP abundance while in the light 
green areas the DM candidate is subdominant. The white regions are excluded since they give a 
too large relic abundance. For very small hidden photon masses annihilation proceeds only via 
the left diagram and is essentially independent of m^' . Therefore the relic abundance is given by 
the kinetic mixing only which itself is determined by the hidden gauge coupling up to a factor k; 
the effect of varying k on the relic abundance is shown in figure [3] where the two extreme cases of 
its range are used. The coloured lines with named labels represent future searches, which as can 
be seen from the plot will probe portions of the interesting parameter space. 

This DM particle can also scatter elastically on nuclei. As it is a Dirac fermion, this process is 
spin-independent and the corresponding cross sections can be compared to the positive observa- 
tions of DAMA and CoGeNT. The results in figurejsjare given for the Standard Halo Model (SMH) 
where part of the CoGeNT allowed region is not excluded by the other experiments (CDMS and 
XENON). The dark matter mass of 6 GeV used in the two top plots is not constrained by CDMS 
or XENON. The band in purple corresponds to the 90% (lighter purple) and 99% (darker purple) 
contours where the correct cross section for CoGeNT can be obtained. At the place where this 
band overlaps with the dark green region, the DM candidate that fits CoGeNT is also providing 
all of the dark matter in the universe. In the larger part however, where the purple band is on 
top of the light green area, the DM particle explains the CoGeNT signal while only contributing 
subdominantly to the total DM. 

For the 7 GeV dark matter the scattering cross section is constrained by CDMS (still below 
the reach of XENON) which is shown as a turquoise line excluding all the parameter space above 
it (where there are dotted vertical lines). Where the CDMS exclusion bound enters the WMAP 
allowed, coloured, region the limit is rescaled as described above to correspond to the appropriate 
dark matter density. However, outside of this region it is not rescaled - the straight line shown 

^This is obviously based on the reasonable assumption that the local DM has the same content of different DM 
contributions as averaged over the whole universe. 
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fjiy, [GeV] 

Figure 3: Dark matter (DM) relic abundance and direct detection (DD) cross section in agreement 
with CoGeNT for a dark matter mass of 6 GeV (top) and 7 GeV (bottom) and the two extreme 
values of the range in which the free parameter k can vary (only for 6 GeV). The Standard Halo 
Model (SMH) has been used. The grey regions are excluded, the coloured lines correspond to 
sensitivities of future planed fixed target experiments. Dark green regions give the correct relic 
abundance, light green is subdominant DM, white is overabundant and therefore excluded. In 
the purple area, the cross sections are such that they can explain the CoGeNT observation. The 
turquoise line is the CDMS bound which excludes all the dashed shaded area above the line. 



corresponds to the behaviour for a constant dark matter density equal to that observed. This 
accounts for the sudden change in gradient. Note that relation (2.4) has a significant effect upon 
the behaviour of this bound. Outside of the WMAP allowed region, i.e. when we are applying the 
CDMS bound for a fixed dark matter density, the corresponding cross-section follows a contour of 
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my, [GcV] 



Figure 4: Scan over hidden photon mass, kinetic mixing and hidden Dirac fermion mass over the 
range from 0.8 GeV to 20 GeV, using the Standard Halo Model. Dark coloured regions indicate 
that the correct relic abundance can be found, lighter colours indicate that the hidden fermion 
is a subdominant dark matter candidate. Green regions are thus simply WMAP allowed, but 
also shown are regions where the direct detection cross section can explain the signal from either 
CoGeNT (purple), DAMA (red) or both at the same time (blue). The parameter k, has been fixed 
to 1. For subdominant DM the scattering cross sections have been rescaled. All points shown in 
the figure are in agreement with all direct detection constraints. The constraint from electroweak 
precision tests is shown as the almost horizontal blue long dashed line, the (model dependent) 
BaBar limit is shown as a dark dashed line and from the muon g — 2 dash-dotted. 



X cx m^/, rather than x m-y which we would find if we were instead keeping gh constant. This 
arises since the cross-section behaves as 

o-DD oc — J- cx (3.5) 

This explains the straight line portion of the CDMS bound in the log — log plots of figure [3| Note 
how this changes when we take rescaling into account: since the thermal-averaged ip-ip annihila- 
tion cross-section multiplied by speed (ciAnn^') for fixed dark matter and hidden photon mass is 
proportional to either gf^x'^ or gf^, which according to equation (2.4) translates into {(TAnnv) oc x^; 
and as the relic density is propotional to l/(a"Ann't^)) we find 

X^ 1 1 
a^(x^- r(x^. (3.6) 
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Hence the rescaled CDMS exclusion bound is approximately a vertical line on the plot, as can be 
seen in figure [3] where the turquoise line meets the green band. 

3.3.2 Scanning over dark matter masses 

To fully examine the parameter space of these models, we performed a scan over the mass of 
the particle ip for values between 0.8 and 30 GeV (in steps of 0.2 GeV) while also varying the 
kinetic mixing and hidden photon mass. The resulting scatter plots are shown in figure |4] for 
three different halo models (for details of the halo models see [66]). For these plots the parameter 
K was fixed to its central value of one. The colouring is as follows: dark shades correspond to 
the dark matter candidate producing the observed relic abundance, lighter shades indicate it is 
subdominant; green regions do not correspond to an experimental signal but are not excluded; 
purple corresponds to explaining the CoGeNT signal; the DAMA signal is explained in the red 
regions; CoGeNT and DAMA signals are explained simultaneously in the blue regions. For the 
Standard Halo Model (SMH), the CoGeNT and DAMA signal regions do not overlap, and this 
is reflected in the absence of blue in our SMH plots. However for the Einasto and NFW halo 
models there is a small region of overlap in mass-cross-section space for the signal regions, which 
translates into blue regions of our plots with those halo models. However, in the following section 
we shall use mostly the standard halo model; the choice of halo has a more dramatic effect on the 
presence (or otherwise) of overlap of the signal region than the allowed parameter space of our 
models. 



4 Analysis of a supersymmetric dark sector 

In this section we describe the results of a scan over the parameter space of the model of section [2| 
constrained by dark matter abundance and direct detection cross-sections. This was achieved by 



implementation of our models in micrOMEGAs ^-50 using LanHEP 72 -76 . We included all of 



the interactions between the hidden and visible sector including the neutralino mixing (described 



in appendix section B.2) and Higgs portal term (described in appendix |E]) which we believe to be 
novel results; as a result, there is some dependence on the visible sector spectrum and couplings. 
Since we were investigating the effects of gravity mediation, and for minimality, we chose the 
visible sector to consist of the MSSM with a Higgs mass above the LEP bound and the lightest 
visible sector neutralino in the range 100 to 200 GeV; the effect of changing the spectrum within 
these ranges leads to quantitative changes of a few percent, but not qualitative ones. 

4.1 Radiative breaking domination 

Here we perform a scan over A5, x ^^^^ ^^7' in order to find parameter combinations which give 
a light dark matter (DM) candidate (target DM masses in the range between 0.8 and 20 GeV) 



which as mentioned in section 2.5 we find to be exclusively a Majorana fermion 01. We insist 



that A5 and the hidden gauge coupling inferred from x remain perturbative; this places an upper 



limit upon x via equation (2.4). We are interested in light hidden gauge bosons, so we choose a 
maximum value of my of 40 GeV. The low energy parameters are found by choosing boundary 
conditions at the high energy scale (10^^ GeV) and running down; this ensures that we have bona 
fide consistent models at the low energy scale, rather than choosing the parameters completely ad 
hoc. This search uses the RGE engine from SoftSUSY j77j . 
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Figure 5: Allowed space of models with radiatively induced breaking. Dark green areas allow for 
the correct dark matter relic density; light green for subdominant dark matter. The Standard 
Halo Model (SMH) has been used, and all constraints imposed except for the SIMPLE exclusion 
limit. Left: k = 1, right: 0.1 < k < 10. 
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Figure 6: Allowed space of models with radiatively induced breaking. Dark green areas allow for 
the correct dark matter relic density; light green for subdominant dark matter. The Standard 
Halo Model (SMH) has been used, and all constraints imposed including SIMPLE exclusion 
limit. Left: n = 1, right: 0.1 < k < 10. 



We then input the results of the scan into micrOMEGAs to obtain the corresponding relic 
abundance and scattering cross sections. Results for kinetic mixing against hidden photon mass 
are shown in figures [5] and |6] for k = 1 and for a scan over k in the range 0.1 to 10. Depending 
on whether the relic abundance corresponds to the total DM abundance the points are shown in 
dark green or in light green if it is subdominant. Clearly allowing for variation in k has a large 
effect on the allowed space of models, in stark contrast to the toy model in section 3.3 where the 
parameter space could be filled without varying k. 

The spin-dependent and spin-independent direct detection cross-sections are shown in figure 
[7| where k has been scanned over and we again use the rescaling procedure for subdominant dark 
matter. The dark matter particle considered in this subsection oi is a Majorana fermion and 
therefore has greatly suppressed spin-independent scattering on nuclei (so there is little chance 
of explaining the DAMA or CoGeNT signals via spin- independent scattering with such a model); 
however, the obtained spin-dependent cross sections are quite large and some are even already 
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Figure 7: Direct detection cross-sections for radiatively induced breaking, using SMH, scanning 
over 0.1 < K < 10. Top: spin-dependent scattering (o"^^) with experimental exclusion contours, 
SIMPLE shown as the lowest-lying, brown, curve; left: spin-independent scattering (o"^^) without 
SIMPLE limit; right: spin- independent scattering (o"^^) with SIMPLE limit. 



excluded by current experiments (the experiments were mentioned in section [s]). The most strin- 
gent constraint arises from the SIMPLE experiment which cuts out many parameter points for 
dark matter masses above ~ 6 GeV; the effect of this is illustrated by showing the parameter scan 
before SIMPLE is included in figure [5] and afterwards in figure [6| Since the spin-dependent and 
spin-independent cross-sections are related, the SIMPLE limit removes a large portion of the pa- 
rameter space with larger values of spin-independent scattering direct-detection cross-sections; this 
is illustrated in figure [7j which shows on the top plot all the found spin-dependent cross-sections 
and the experimental bounds, then on the left and right plots the corresponding spin-independent 
cross-sections without and with the SIMPLE bound respectively - the missing models on the right 
plot are those that lie above the SIMPLE exclusion contour for spin-dependent scattering. 

As mentioned above, the fact that the hidden sector dark matter candidate is a Majorana 
fermion leads to extremely small spin- independent scattering cross sections. They do however 
obtain a contribution from the Higgs portal term, which in supersymmetric theories is always 
present. We describe this in detail in appendix [Ej where we also derive a simple approximation 
for the contribution of the Higgs portal term which agrees well with the results seen in figure [7| 

~10-«om= X (-^)\l^y- (^)'. (4.1) 

\mj\f -\- tUq-^^ / vu.UUi/ \ my J 

There is also a somewhat smaller and more spectrum-dependent contribution from squark ex- 
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change. Here the spin-independent nuclear cross-sections are very similar for protons and neu- 
trons; that for protons is shown for direct comparison with the next subsection. 

Our results in this subsection are largely independent of the halo model applied to the spin- 
independent scattering limits, as the corresponding cross-sections are much below the experimental 
reach. 

4.2 Visible sector induced breaking 
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Figure 8: Allowed space of models with visible-sector induced breaking, showing dark green areas 
where the correct dark matter relic abundance can be found and light purple where the CoGeNT 
signal can be explained with a subdominant dark matter candidate. The Standard Halo Model 
(SMH) has been used, and all constraints imposed including the SIMPLE exclusion limit. Left: 
K = 1, right: 0.1 < k < 10. 



Here we implement a scan for visible sector induced breaking, by scanning over parameters at 
the low energy scale. As in the previous subsection, we insist on perturbativity for Xs and gh, 
and take a maximum value of my of 40 GeV. However, the soft supersymmetry breaking masses 
are chosen to be small, as may be induced in gauge mediation or sequestering of the hidden 
sector. Phenomenologically then, the results of this subsection can be considered to be a detailed 
examination of the model of |T9], but with a large gravitino mass and kinetic mixing respecting 
the relation (2.4). 



As mentioned in section 2.5, we can have either of two dark matter candidates, depending 
on the particular low energy parameters: either the Majorana fermion 6i or a Dirac fermion 07. 
For both cases we again use micrOMEGAs to compute the relic abundance and the scattering 
cross sections. The space of models in the kinetic mixing-hidden photon mass plane (all points 
shown are in agreement with all direct detection exclusions, including SIMPLE) are shown in 
figures [8] and [9} the colour code is identical to the scatter plots in the toy model of section |3.3| 
and the different experimental constraints are explained in the constraints section |3] (sensitivities 
correspond to DAMA and CoGeNT 90% and 99% contours as mentioned earlier). Figure [s] 
demonstrates the expansion in the parameter space by allowing a varation in k; both k = 1 
and 0.1 < K < 10 are plotted, for the Standard Halo Model (SMH). Figure [9] shows the effect 
of changing the halo model: here we find subdominant explanations for separately DAMA and 
CoGeNT signals, and both simultaneously, which appear as light red, purple and blue regions 
respectively; simultaneously explaining both signals is only possible for certain halo models other 
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than SMH. 
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Figure 9: Allowed space of models with visible-sector induced breaking, scanned with 0.1 < k < 10, 
showing different halo models including the SIMPLE limit: top: NFW halo model, left: Isothermal 
halo model, right: Einasto halo model. Here the red region shows the space explaining the DAMA 
signal, the purple region explains the CoGeNT signal, and the blue explains both DAMA and 
CoGeNT, all signal regions having a subdominant dark matter density. 



The resulting parameter points in figures [8] and [9] show a very similar behavior to the toy 
model (see figure |4]) as here the dark matter candidate can also be a Dirac fermion. The main 
difference is that the models here never permit annihiliation via the t-channel diagram - since the 
dark matter particle can never be heavier than the hidden gauge boson. Therefore the lower part 
of the plots is, in contrast to the toy model, empty as there it was filled by dark green points 
finding the correct relic abundance lying either just above the threshold for t-channel annihilation 
or the s-channel resonance (as we scan over the dark matter mass these resonances move through 
the plot through different values of m^i and the whole range is covered). The coarser grid and 
small holes in the current scatter plots compared to the toy model arise from the fact that the 
parameter space can not be scanned as continuously as for the toy-model. 



The spin-dependent cross-sections are shown in figure 10 , which are only appreciable when the 



Majorana fermion oi is the dark matter candidate. The spin-independent cross-sections for the 



Standard Halo Model are shown in figure 11, where the effect of scanning over the parameter k 
on the parameter space is illustrated by showing plots with k = 1 and 0.1 < /t < 10. There are 
two disjoint regions corresponding to whether the dark matter candidate is the Majorana fermion 
oi (lower region) or Dirac, 07 (upper region). As in the radiatatively-induced breaking case the 
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Majorana fermion has a small spin-independent cross- section of 10"^^ to 10"^^ cm . In contrast 
to this, the spin-independent scattering cross section of the Dirac fermion 07 is in the range of 
current direct detection experiments, and may explain the signal due to CoGeNT and DAM A via 
a subdominant dark matter component. Hence we show the effect of the different halo models in 



figure 12, where simultaneous explanations of both signals are possible. 
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Figure 10: Spin-dependent cross-sections for visible-sector induced breaking, with the exclusion 
contours from direct detection experiments shown (the lowest lying, brown, line is due to SIMPLE). 
Here only the Majorana dark matter candidate oi is shown (the cross-sections due to an o^ 
candidate are too small to appear). The SMH has been used. Left: k = 1, right: 0.1 < k < 10. 
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Figure 11: Spin-independent cross-sections for visible sector induced breaking using the SMH 
and including SIMPLE limit. The signal contours from CoGeNT (purple) and DAM A (Red) 
are shown; the exclusion limits from CDMS and XENONIOO are shown as the coloured curves. 
The plot splits into two disjoint green areas: in the upper, Dirac fermion o^ is the dark matter 
candidate, while in the lower one it is the Majorana fermion oi. Left: k = 1, right: 0.1 < k < 10. 



The Dirac fermion 07 has almost no spin-dependent scattering on nuclei, so the effect of the 
SIMPLE exclusion limit only affects the Majorana fermion (lower) regions in the spin-independent 
scattering plots 11 and 12 while the parameter regions that are interesting for spin-independent 
scattering experiments and can explain the DAMA and CoGeNT signals remain untouched. The 
spin-independent cross-sections are those on protons; as described above, the interaction through 
kinetic mixing couples almost exclusively to the charge of the nucleon, so the spin-independent 
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scattering on neutrons essentially vanishes (for the Dirac fermion 07). For the case of the Majorana 
fermion oi , spin- independent scattering on protons and neutrons is of roughly equal magnitude as 
it preceeds via the Higgs portal and squark exchange (as in the previous subsection), but the plots 
also show only the cross section on protons. Spin-dependent scattering is also shown in figure 10 
on protons only but as in the previous subsection the cross sections on neutrons are almost the 
same (for the Majorana case). For the Dirac fermion dark matter the spin-dependent scattering 
essentially vanishes both for protons and neutrons. 
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Figure 12: Spin-independent cross-sections for visible sector induced breaking including SIMPLE 
limit, scanning over 0.1 < k < 10, and showing three different halo models. Top: NFW, left: 



Isothermal, right: Einasto. The signal and exclusion regions are shown as in figure 11; as there, 
the upper and lower green regions in each plot corresond to dark matter candidates 07 and oi 
respectively. 



5 Conclusions and outlook 

We have presented what we believe to be the first detailed examination of the dark matter relic 
abundance and direct detection cross-sections of a complete supersymmetric dark force model, 
emphasising the natural supersymmetric relationship between kinetic mixing and the hidden gauge 
coupling. In particular, we have included running from high energy gravity-mediated boundary 
conditions, and shown that interesting and viable models exist, in contrast to prior expectations. 
We have also examined the effect of neutralino mixing and the Higgs portal term, showing that 
the latter can contribute a small spin-independent cross-section for Majorana fermion dark matter 
candidates. We examined the model in the cases of both radiative and visible sector induced hidden 
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gauge symmetry breaking, and demonstrated the stark phenomenological constrasts between the 
two. 

While the model can be used to explain the current dark matter signals observed by DAMA, 
CoGeNT and CRESST, this is not plausible in the case of radiative-induced breaking relevant for 
gravity mediation, where our motivation was to show that simple dark sectors are not excluded 
- the hidden U(l) may instead be detected in fixed target experiments, particularly if the hidden 
photon cannot decay to hidden matter (as in the reasonably generic case when the dark matter 
particle has mass near that of the hidden photon). However, this is certainly plausible if the model 
were to be extended, for example by allowing a supersymmetric mass for the singlet. 

Wc hope that this work has paved the way for more detailed analysis of other supersymmetric 
dark sectors. In addition, there are several further possible avenues of work within the current 
model, for example by including the full loop corrections to the effective potential, and constraints 
from indirect dark matter searches (which is a work in progress although we believe them to be 
much less stringent than the direct searches). It would also be interesting to compare the signal 
from CRESST with those from DAMA and CoGeNT in the context of these models. 
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A Renormalisation group equations 



Here we present the two-loop renormalisation group equations for the hidden-sector parameters 
as = A|/47r, = gf^/An, Mx, m\, m^,mg, As, neglecting the effect of kinetic mixing. We define 
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t = log//: 



^ =^[2as{3as - 4ah)] + -^[4as{-Sal + 2asah + 8a|)] 

^ =^ [6«5^5 + Sa^MA] + [8 4 ^rn,. _ 

1 



4aftMA] + — -2 [8a/,MA(4< - as) + Sa^Qs^s] 



+ 



^as{ah - a|)(m| + ml + ml + 2^|) + 8afta5(2M| - 2MaA5 - A|) 



(47r)2 

=^ - [2as{m's + m\ + w?L + A|) - 8M|a,j ± 2ah{m\ - m^)] 
1 



dm^ 1 



+ 



(47r)2 _ 



+ 8Q!^(m+ + rrii) ± 4a(2Q;/i - Q;5)(m+ - m?.) 



(A.l) 



B Spectrum of the model 



In this appendix we present the details of the low energy features of the model W = XsSH-i-H-. 
Once supersymmetry and the R-symmetry is broken the potential is 

V =|A5p(|5i?+p + \SH_\'^ + \H+H_f) 

+ §{\H^\'-\H-f-0' 

+ ml\H+\^ + ml\H_\^ + ml\S\'^ 

+ {XsAsSH+H_ + c.c.) (B.l) 

We are assuming that no /j,, terms are generated; these would introduce new scales into the 
theory. Although they could conceivably be generated by a Giudice-Masiero mechanism in analogy 
with the visible sector, we shall neglect this possibility. 

B.l Scaleirs 

Defining A = ^A|^ — rn^X^/gf^, we have mass matrices in the basis of 

which implies masses for the two components of v^S'aA/A^, at this level. The "massless" mode 
is the Goldstone boson that becomes the longitudinal component of the massive gauge field. The 
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{H-,hI,S,S^) system is more complicated; we find a mass matrix of 



hnl i7_ s) 



( + + m?. 
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(B.3) 



The above theory for non-zero ^ has a minimum at = A/A5 provided that A is real, and 

<m^ + + m| + 2^^ (B.4) 
<{mL + + A2)(m| + A^) - l^^pA^. 

In the case As = 0,m\ = mi < this translates simply to the condition that A| > 2g^. 



B.2 Fermions and neutralino mixing 

The fermion mass matrix in the basis (A, h-\-,h-, s) (neglecting the kinetic mixing of the gaugino 
with the neutralino) is given by 



Mf 



( Mx 









my 
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(B.5) 



However, to properly compute the dark matter density we should take mixing of the fermions 
with the neutralino into account. The fields h-,s form a Dirac fermion that does not mix with 
any other fields. There will however be kinetic mixing of the bino with the hidden gaugino, and 
possibly mass mixing; writing these fields before the mixing as respectively b, A we can define 




A =A - teb, (B.6) 

which then allows us to write the full neutralino mass matrix in the basis {Bq, Wq, h^, h^, A, 
including the standard Majorana masses for Bino and Wino Mi^2 but also a potential explicit 
mass mixing term jO D —mxbX as 



/ ^jMi-SeUix + s'^.Mx] -Mzswc^/ce MzswS/3/ce ^^[mx-SeMx] -m'^t,\ 
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-MzSwCi3/Ce 
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MzSwSp/Ce 


-MzcwS0 


-A 











^^[mx-SeMx] 











Mx 

























(B.7) 
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C Kinetic and mass mixing 



Here we review the diagonalisation of the gauge fields and the subsequent couphng of the physical 

gauge bosons to matter fields. 
Consider the lagrangian 



0" 



1 



+ T—^T,Xf.X^ + Iv^gyB^ + ^^ X^ - g2W^f 



Then we rotate 



(C.l) 



B^ ^cwA^ - sw{c^Z^ + s^7^) + 



0^ 



:(-S<^Z^+C07;[,) 



(C.2) 



where cyi/,syY are the usual cosine and sine of the weak mixing angle respectively, C(i),S(i) are the 
cosine and sine of an an angle (j) to be determined below so that 



1 



1 



3w 3b 



+ Z^ 



+ 7 



XS0 



92Cwc^Jw - \swC4> + )9yJb 



0^ 



3x 



92Cws^jw + ( 



C0X 



SWS,j>)gYjB + 



9xc^ 



0^ 



(C.3) 



We find, defining x = m? /M"^: 



tan 2^ = 5 5 — ^ 



sin ffl = h 

1 — X 



(C.4) 



In terms of these, we have the eigenvalues where m+ corresponds to the physical Z mass, m_ 
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to the physical hidden photon mass: 

1 



9 Ml 
cos^ a 



cos^ a 



mi =Mi 



1 ~ 1 
m_ =m 



1 — X 



1 — X 



+ 



(C.5) 



Thus the masses are only shifted at order x^- 
Note that we can also write 



m 



rrij 



~ 2 

m 



m_ =Mz{s^ + swteC(j,) + ^c; 



(C.6) 



and, defining x = m?_/m\ ~ x we have 



tan ( 



-(1-X)± J(l-x)2-4sy2^ 



2swtfX 



(C.7) 



D Goldstone boson mixing 



Here we consider the eflFect of the mixing on the Goldstone bosons that are eaten, and what 
happens to the other fields. Assuming that the visible sector is the MSSM, and taking the hidden 
sector to be broken by a single complex scalar C, we have 



L D|9^C - ighX^C\' + \d^H^^ - -{gyB^ - ^sO^^P + + ^{gyB^ - g2W''^)H^^\' + Vhid{C) + V,is{H) 



0\ ij0|2 



-l{d^Ci-mhX^f + \{d,CRf + \ 



d^ic^Hs - - -{gyB^ - g2Wll)vs^ 



+ 



d^{spH^ + CfsG'') + ligyB^ - g2W?,)vcp 



+ ... 



Z)-{d^Ci-mhX^f + 



ev 



S2W 



+ ... 



(D.l) 



Clearly the masses of the neutral Higgs (both visible and hidden) are unaffected by the mixing, 

(D.2) 



but the pseudoscalar will be. Using (C.2) and 



CWW^I - SwB^ =(C0 - SwteS^)Zfj_ + {s^ + SwteC^h'fj. 



we obtain 



(D.3) 
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Thus 



m+Gz =Mzoic^ - swteS(j>)G^ - ruh — Ci 
m-Gy =Mzo{s^ + swteCcj))G'^ + vrih—Ci 



and 



Ci = 



Mzo 

Ce 

rrih 



mhC^m+ mhS^m- 
(j-Z H t^V 

C(f,m+Gz + s^iJi-Gy 

- (s</) + swteC(j>)m+Gz + (c0 - swteS(j,)m-Gy 



(D.4) 



Mzorrih [ c, 
1 



We can also write this as 



G° = cos ipGz + sin V^G-y/ 
Cj = — sin V'Gz + cos '^Gy 



tan^ =(tan^) 



E Higgs portal mixing 

With the Lagrangian density 



X 



we need the D-term mixing; we write Dy = —9y 4''''^4': B>x = —Qh J2q 
in the absence of mixing, then we have 

>C 3^^y + - xDyDx - DyDy - DxDx 



which leads to 



Dx 



1 



;{Dx + xDy) 



DY=Y^iDY + xDx) 



and thus 



Portal = 



2 



1)2. + + 2xDxDy 



DxDy. 



(D.5) 



(D.6) 



(E.l) 

x(p &s the D-terms 
(E.2) 



(E.3) 



(E.4) 
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The relevant part of the potential for us involves the Higgses; we can write the portal term as 

X „ _ /irr |2 ITT |2\ /-^ I rr |2 I rr |2 



^Portal 



1-x 
xgygh 



29Ygh{\H+f 



\H-f){l\H,f-l\H,n 



(E.5) 



i\H+f-\H^\') 



2' ^2 



Immediately we can extract the effective Fayet-IIiopoulos term: 

Vpj=^j{\H^f-\H^f-0' 



2n 



^^{xl9h){DY) 
={x/9h)9Y^ cos 2/3. 



(E.6) 



However, we can also extract the Higgs mass mixing. Writing 



1 

1 



[s^v + hl + i{cpA - spG° 



V2' 



cos a — sm a 
sin a cos a 



H 
h 



(E.7) 



and using the standard shorthand = cos /3, Cq = cos a, c^+p = cos(a + /?) etc, we have 
F D - — c/yc/fel^ + ^(xR + zx/)! 



X 



0mi2 



D 9Y9h7;v+vxR 



=—9Y9hV 



V2A 



= — 2M^s^ymy — 2xsvF-^z'T^7'- 
So then we must redefine our Higgses: the mass mixing matrix in the basis h,H) is 



(E.e 



•^Higgs - 



-s^+pMl ml 



(E.9) 
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To first order in M'^/m\^ tliis is diagonalised via 
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ml-ml 
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(E.IO) 



Since we often find rnjj ^ m| however the above will usually reduce to mixing between the 
hidden and lightest Higgs. In this case we can approximate 
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1 u 
-u 1 
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Mzm^ 
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(E.ll) 



Thus very roughly u ^ x^wm-y' / nih for large tan^S. 
E.l Spin-independent nucleon cross-sections 

Here we would like to estimate the cross-sections for our Majorana hidden fermion Oi 



oi 
oi 



on nucleons that take place via the Higgs portal term. The Higgs portal leads to an effective 
four-point interaction 



CD fN{OidiNN] 



(E.12) 



and we can consider different fp, fn for protons and neutrons respectively. Consider that the dark 
matter particle is a Majorana combination of A and h+ fermions, oi « cos 9ih+ + sin^iA. The 
coupling to the Hidden higgs is via the kinetic vertex 



C D - V2ghHl{h+~^) + c.c. 

D - Qh cos 01 sin 6ixr[{oiOi) + (oioi)] 

D — gh cos 9i sin OixrOiOi 



(E.13) 



oi 



where now Oi is in Dirac form, Oi = - . Now let us write the coupling of the MSSM Higgs 
_ \oi J 

to nucleons as —aNhNN. 

L D -5/1 cos 01 sin0i[x'^ + 'u/i']Oi6i - a^\-ux'ji^\{\NN - ^ml{x'iif - ^ml{h'f 
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(E.14) 
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So then our direct detection amplitude should be approximately given by the coefficient above. 

The coupling of the MSSM Higgs to nucleons is determined by its coupling to quarks. These 
come from 



£ D - YuHuqu - YoHDqd 

D —hijnijuu) —himodd). 



(E.15) 



Then we use that 78 



{N\m,qq\N) =m„4^^ 



{N\mQQQ\N)Q=c,b,t =^^N 

2 

= 27"^^ 



q=u,d,s 



1 _ i?{JV) 



(E.16) 



to give 



Cq, 

Oat = m^r 



H niN- 
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If we consider large tan /3 with /3 = 7r/2 — 5, q ~ —6 we have Sc 



-S,Ca ^ l,Ci3 ^ 5,Si3 ^ 1 



and so 
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[^-^) + Jtu - JTd - J: 



Ts 



Let us then simply define 



OAT 



rriN 



(E.18) 



There are large uncertainties in the value of /'•^•*, however it is approximately equal for protons 
and neutrons, and varies from about 0.03 to 0.44. We can then take an approximate value for the 
amplitude for a Majorana fermion scattering on nucleons to be 



/at PS sm26'iSa+/3XSH/ — o — k J ■ 



(E.19) 



Taking the large tan/3 values, ruh = 115 GeV and /'•^•' ~ 0.1 we obtain 



/tv « 3 X lO-y(GeV)-^ x sin2e 



Xsw 
0.001 



(— ) 

\ my J 



(E.20) 



This is clearly a very small effect. This translates into a cross-section for scattering on a single 
nucleon of 



_SI, Portal 
N 



7r{mN + rrioi 



il X 10""^^ cm^ X 



ol 



sin2(9i 



0.001/ 



/GeVy 
\my ) 



(E.21) 



This corresponds well to the values that we find in the plots. 
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The above will be supplemented by contributions from s-channel squark exchange. Very 
roughly for these we have effective four-point interactions of the Majorana fermion with quarks 
with coupling fg ~ '£[t\Ui- where b = J2n ^bo^On is the Bino, which mixes most strongly with 



the lightest hidden state. By considering the mass mixing matrix in section B.2 we can conclude 
that (in the absence of direct mass mixing) the mixing is simply ~ X^A; so we have C/^^^ ~ x, 
and thus 



and thus 



SD,Squark ^^-48 2 "^li \ f SwX ^ ( WOGeV\ 



4 



^i.,.quarK ^^g-^^Cm^ X , """^ ' . (E.23) 
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